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1. Motivation and Objective 

The primary objective of the Center for Modelling of Turbulence and Transition 
(CMOTT) is to further the understanding of turbulence theory for engineering 
applications. One important foundation for this research is the establishment of 
a data base encompassing the multitude of existing models as well as newly proposed 
ideas. The research effort described in the next few pages is a precursor to an extended 
survey of two equation turbulence models in the presence of a separated shear layer. 

Recently, several authors' have examined the performance of two equation models 
in the context of the backward facing step flow. Conflicting results, however demand 
that further attention is necessary to properly understand the behavior and limita- 
tions of this popular technique, especially the low Reynolds number formulations. 
The objective of this research is to validate an incompressible Navier Stokes co e 
for use as a numerical test-bed. In turn, this code will be used for analyzing the 
performance of several two equation models. 

2. Work Accomplished 

To date, the validation of the incompressible code DTNS is complete. The details 
of this validation study are documented in reference^]. The code is based upon 
the pseudo-compressibility technique and incorporates the approximate factorization 
scheme for time integration. Two laminar benchmark flows are used to measure the 
performance and implementation of the numerical methods. The classic solution of 
the Blasius boundary layer is used for validating the flat plate flow, while experimental 
data is incorporated in the validation of backward facing step flow. 

An initial result for the standard high Reynolds number k-e equations has also been 
calculated to demonstrate an initial performance level of the solution technique for 
the turbulence equations. 

2.1 Numerical Method 

The absence of a variable density in the continuity equation governing incompress- 
ible flow complicates the numerical integration procedure. One solution technique 
that has been well received is that of pseudo-compressibility. This idea was first 
put forward by Chorin[2] and enables the equations to be solved using the prim- 
itive variables. Recently, Chang and I<wak[3], Rizzi and Eriksson[4], Kwak and 
Chakravarthy [5] , Michelassi and Shih[6], and Turkel[7] have found this method suit- 
able for resolving incompressible flow. This particular implementation has been vali- 
dated by Gorski[S-10] for several different benchmark flows. 
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Examine the system of equations involved in the pseudo compressibility method 
and notice that they differ from the incompressible Navier Stokes equations by the 
addition of a time dependant pressure term in the continuity equation. 
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Here, x and y are the independent variables and Re refers to the Reynolds number. 
The constant ft is known as the pseudo compressibility parameter. This system 
is hyperbolic in nature whi le the incompressible flow equations are elliptic. The 
pseudo sound speed, c = a/u 2 + /5, is governed by the value of the parameter /3, 
whereas the physical sound speed is infinite. Chang and Kwak [3] have shown that 
for ft > 0 the finite speed pseudo waves vanish as time progresses and yield the proper 
incompressible solution at the steady state limit. It is through this parameter P that 
the convective and acoustic waves are decoupled, and thus convergence is governed. In 
choosing an optimum value for this parameter, the goal is to avoid giving the viscous 
effects time to react to the passing of the nonphysical transient pressure waves. Thus 
a lower bound on the acoustic speeds translates into a lower bound on fi. However, 
an upper bound on ft is strictly scheme dependent. 

The approximate factorization is rather straight forward. Equations (2) can be 
rewritten as 
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or in generalized coordinates U(x,y),y(x,y)) as seen here in equation (3). 
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The third order accurate Chakravarthy and Osher TVD scheme is used to discretize 
the convective terms while central differencing is used for the viscous ifluxes^. T 
the approximate factorization of this method is expressed in terms of flux Jacobians 

A and B as 
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where V is the cell volume, and At = » the timestep. The matrices A, , and IB. 


wnere v is me ecu vuiu.uc, — - 14 VJ . rQl 

are the viscous flux Jacobians. The interested reader is encouraged to see reference^ 

for more details. . 

The k-e equations are solved in a similar manner as that seen m equation 2. 
Consider the standard high Reynolds number form of the transport equations of the 
turbulent kinetic energy and dissipation of turbulent genetic energy as in equation 

(5): 
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and the production of turbulent genetic energy defined as shown below. 

P = vt [2(4 + v y) + ( u y + v *) 2 ] 

Again, we can translate these equations in a manner similar to that of equation (3) 
and create the following system for generalized coordinates: 
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This can be solved via approximate factorization as shown here in equation 9. 
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Again, the scheme is written in terms of the turbulence flux Jacobians and the 
approximate source Jacobian, H 4 . The coefficients and a-q are provided to enhance 
the diagonal dominance of the implicit method. While the convective flux terms are 
linear, the linearization of the source terms is not straight forward. Furthermore, the 
low Reynolds number formulation of equation 5 includes other nonlinear correction 
terms in the source vector S 4 . The interested reader is encouraged to see the recent 
paper by Michelassi and Shih[6]. 


2.2 Discussion 

Overall the performance for laminar flow documented in reference[l] is very en- 
couraging. The flat plate results shown agree with the theoretical Blasius solution 
for several different grid configurations. The laminar back facing step flow was also 
well resolved both in terms of the primary recirculation zone reattachment length, 
and corresponding velocity profiles as seen in figures 1 and 2. However, performance 
of the two dimensional flow code drops off dramatically when the Reynolds number 
corresponding to the onset of 3D flow is exceeded. Armaly et al [11] observed this 
phenomenon occur at a Reynolds number of 400, based upon the hydraulic diameter 
of the inlet channel. Again, the numerical results bear this finding out. 


Figure 1 Reattach ment length predictions — comparison with experiment [11]. xj, X2, and X3 are the 
primary reattachment, secondary detachment, and secondary reattachment lengths, respectively. 
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Figure 2 ReD=389 velocity profiles — comparison with experiment [11]. 
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The current thrust of effort is in establishing a comparison between the standard 
high Reynolds number form k-e model and a version which incorporates the low 
Reynolds number corrections, such as that of Jones and Launder [1 ]. is mo e 
is particularly convenient in that the damping functions are independent of the 
parameter y + Thus, multiple walls bounding the computational domain do not need 
special attention. An initial calculation using the high Reynolds number approach 
agrees with the published results of Speziale and Thangam[13]. However, the Jones 
and Launder model is proving to be quite challenging. 


3. Future Plans 

Eventually, several two equation models will be examined for their behavior in the 
back facing step flow. The following is a partial list of the two equation models of 

primary interest: 

1. Jones and Launder 

2. Chien 

3. Yang and Shih 

4. Shih and Lumley 

The current research effort will result in an increased understanding of the need 
for proper modelling of the near wall region. However, the increase in accuracy 
maybe offset by an unwarranted increase in effort. Further numerical experiments 
are necessary before any conclusions can be drawn. 
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